Semi-Hard Scattering Unraveled from Collective Dynamics by Two-Pion Azimuthal 
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Elliptic flow and two-particle azimuthal correlations of charged hadrons and high-pr pions 
(j>T > 1 GeV/c) have been measured close to mid-rapidity in 158A GeV/c Pb+Au collisions by the 
CERES experiment. Elliptic flow (V2) rises linearly with pr to a value of about 10% at 2 GeV/c. 
Beyond pr ~ 1.5 GeV/c, the slope decreases considerably, possibly indicating a saturation of V2 at 
high pt- Two-pion azimuthal anisotropics for pr > 1.2 GeV/c exceed the elliptic flow values by 
about 60% in mid-central collisions. These non-flow contributions are attributed to near-side and 
back-to-back jet-like correlations, the latter exhibiting centrality dependent broadening. 
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Late stages of ultra-relativistic heavy-ion collisions are 
characterized by strong collective transverse expansion. 
An important signature of collective dynamics in non- 
central collisions is elliptic flow manifesting itself in an 
azimuthal anisotropy of particle yields with respect to the 
reaction plane yj, l2( . It is driven by anisotropic pressure 
gradients built up during the early stage of the collision 
in the geometrically anisotropic overlap zone Q. The 
degree of equilibration achieved during the subsequent 
evolution might in principle be gauged by comparison to 
hydrodynamic models which indicate an upper limit of 
elliptic flow. Its magnitude depends on initial conditions 
and the equation of state (EoS), but also on system life- 
time |4| . It had been argued |5j that hydrodynamics was 
unable to describe elliptic flow data at SPS energies while 
it accurately describes RHIC data 0. 

Two-particle azimuthal correlations are sensitive to 
collective flow but might reveal, particularly at large 
transverse momentum (pt), also relics of primary scatter- 
ing in the semi-hard sector which are masked in inclusive 
Pt spectra at SPS energies by the Cronin effect Q. It is 
the purpose of this Letter to demonstrate that collective 
flow and semihard scattering can be disentangled by mea- 
surements of azimuthal anisotropics. We present data (i) 
of elliptic flow for charged particles and for identified pi- 
ons up to pt ~ 3 GeV/c, and (ii) of two-pion correlations 
at pt > 1.2 GeV/c. We pursue a statistical analysis 8] 



which conjectures that the observed anisotropies of non- 
flow nature are due to dijet-like correlations. An inter- 
pretation in terms of resonance decays is unlikely in view 
of the high invariant mass (rs 2.5 GeV/c 2 ) required. 

CERES has acceptance close to mid-rapidity 
(2.1 < r\ < 2.65) with full azimuthal coverage. Charged 
particles (/i^) are tracked by a doublet of silicon- 
drift detectors (SDD) before, and a multiwire propor- 
tional chamber behind the magnetic field used for mo- 
mentum determination. Identification of pions with 
p > 4.5 GeV/c is performed by two ring-imaging 
Cherenkov detectors (RICH). Pion momenta are deter- 
mined from the ring radii. We analyzed 43 million 
158 AGeV/c Pb+Au collisions taken in 1996 at the most 
central (26±1.5)% of the geometric cross section, a geo . 
We divide the triggered events into six contiguous cen- 
trality bins labeled by the percentage of a geo of their 
upper edges, 26%, 21%, 17%, 13%, 9%, and 5%, respec- 
tively. The centrality is determined by charged particle 
multiplicity (A^) in the range 2 < rj < 3 covered by 
the SDDs. Numbers of participants, N part , and binary 
collisions, N rn u, were calculated from a nuclear overlap 
model [lf| neglecting fluctuations, resulting in a total 
inelastic cross section of 6.94 barn. 

Conventional elliptic flow analysis is based on the az- 
imuthal particle distribution with respect to the orienta- 
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FIG. 1: Centrality dependence of V2 for (a) h*, 0.5 < 
p T < 2.5 GeV/c, and (b) tv ± , p T > 1.2 GeV/c. Statistical 
errors are within symbols, (c) pr dependence of V2, corrected 
for BEC effects. Here, the centrality corresponds to three left- 
most bins in (a),(b) combined. Hydrodynamical calculations 
with phase transition at T c = 165 MeV are shown for kinetic 
freeze-out temperatures of 120 MeV and 160 MeV. 



tion \I/ of the reconstructed event plane (EP), 
dN 



d(<£-*) 



= A[ l + 2v' 2 cos(2(0- *)) ]. 



(1) 



Measurements of the EP and of the particle anisotropy 
are obtained from different subsets of the same data sam- 
ple: the 4> acceptance is divided into 100 slices and every 
fourth slice is combined into a subsample. We avoid au- 
tocorrelations by using non-adjacent ^-slices for v' 2 and 
EP measurements. Non-uniformities in \P arc removed 
by standard procedures [ill ]. Depending on centrality, 
the r.m.s. of EP resolution is 35-40 degrees. The Fourier 
coefficient v 2 is obtained after correcting for the EP dis- 
persion, v 2 — K v' 2 . The correction factor for sample i, 
JCi = (cos(2(4' 4 — Vf^))) -1 , involves the unknown angle 
of the true reaction plane but can be expressed by 
the measured sample differences ^ — ^ . This spread in 
K, factors dominates the systematic uncertainties in V2- 
Resulting absolute systematic errors vary with centrality 
between 0.5% and 1.5%. 

The flow results are presented in Fig.^ With central- 
ity, v 2 decreases almost linearly as expected from hydro- 
dynamics. The larger v 2 values for tt^ reflect their larger 
(p T ) ~ 1.45 GeV/c compared to 0.70 GeV/c for h ± . The 
Pt dependence of v 2 measured for the first time at the 



SPS up to 3 GeV/c is shown in Fig. ^c). The data are 
averaged over the three most peripheral centrality bins 
and correspond to (13-26)% of a geo ((N c h) — 190). It 
is seen that v 2 rises about linearly with px to a value of 
about 10% at 2 GeV/c. Beyond p T « 1.5 GeV/c, the 
slope decreases considerably, possibly indicating a satu- 
ration of v 2 at high px as observed at RHIC Q . 

The data in Fig. QJc) have been corrected for the ef- 
fects of Bose-Einstein correlations (BEC) 12] with input 
from The corrections vary between -15% of v 2 at 

p T = 0.25 GeV/c and +10% at p T > 1 GeV/c. Since the 
procedure becomes questionable for central collisions, the 
data in Figs. IHa,b) were left uncorrected. We may com- 
pare to recent results of NA49 at (12.5-33.5)%. Af- 
ter correcting for different centrality, these are still about 
15% larger in the range 1.0< pr <1.5 GeV/c. 

Results of hydrodynamical calculations || using an 
EoS with a first order phase transition to quark gluon 
plasma at T = 165 MeV and freeze-out at Tf = 120 MeV 
are considerably above the data. If the hydrodynamical 
evolution is terminated already at Tf — 160 MeV [lij . 
good agreement with the elliptic flow data is reached. It 
should be noted, though, that in this case the px spectra 
of protons are too steep (not shown). Conversely, the 
Tf— 120 MeV calculation reproduces the spectra. Pos- 
sible explanations include incomplete thermalization [lfij 
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FIG. 2: Anisotropics V2 from EP method for pions (circles) 
and yfp from two-pion correlations (triangles), (a) Central- 
ity dependence for full azimuth (open triangles) and a range 
restricted to |A0| > 0.6 rad (closed triangles), (b) pr depen- 
dence of V2 for the top (13-26)% of a geo , {N ch } = 190. 
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and viscous effects [l7j . 

We turn to the measurement of two-pion angular cor- 
relations at pj- > 1.2 GeV/c which are written in terms 
of the relative azimuthal angle Acf> as 



~ 2 
■o 
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dA0 



S [l + 2pcos(2Ac/>)]. 



(2) 



From the data we determine the second Fourier coeffi- 
cient p which for pure flow is equal to v\ 11]. The 
two-pion yield is corrected for single-track reconstruction 
efficiency which is determined by embedding simulated 
tracks into real events. This correction varies between 
16% and 9%, depending on N c h- At small opening an- 
gles, overlapping rings in the RICHes cause a drop in 
pair reconstruction efficiency. To be less sensitive to a 
Monte-Carlo (MC) correction, pairs with track separa- 
tion AO < 20 mrad in polar angle were discarded. This 
cut reduces the pair efficiency loss by a factor of four 
while keeping still 60% of statistics. Corrected opening- 
angle distributions reveal strong a strong anisotropy with 
maxima at A(j) ss and tt (see below). The procedure 
is supported by the fact that anisotropics remain essen- 
tially unchanged after correction whether or not the AO 
cut is applied. 

In Fig. Efa) is shown that the ^fp values from the two- 
particle correlation are systematically larger than the i>2 
coefficients. It is unlikely that the i>2 values are signifi- 
cantly reduced by a possible bias on EP reconstruction 
by the high-p^ particles in view of their very small abun- 
dance («10 -3 of h ). We have analyzed the full az- 
imuth applying the A9 cut and MC correction for pair 
efficiency, and alternatively a range \A<f>\ > 0.6 rad with- 
out these remedies. For the first centrality bin ((21-26)%, 
(N c h) = 147), the value of ^fp exceeds v-i by 70% for full 
range in azimuth; accounting for BEC effects, this ex- 
cess is reduced to about 60%. The anisotropy in the 
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FIG. 3: Two-pion opening angle distributions for 
p T > 1.2 GeV/c (a) for the first centrality bin, (21-26)% 
and (b) for the fourth, (9-13)%. A cut A8 > 20 mrad and 
corrections for close-pair efficiency losses are applied. Full 
line shows Gaussian fits to semi-hard components on top of 
flow- modulated background (dashed line). 



O Af:0 
1.5 -• A ^ = 7t 



0.5 



DQ O O o 



250 500 750 1000 

N„„„ 




250 500 750 1000 

N„„„ 



FIG. 4: Centrality dependence of the Gaussian widths of the 
correlation peaks at A<j> = 0, n (a) and of the areas under 
the Gaussian peaks, both from fits (b). The most central 
points for A(j> = it combine the fourth and fifth centrality 
bin, (5-13)%. The loss in pair acceptance due to the cut 
A6 > 20 mrad has not been corrected for. 



restricted A<j) range is similar in magnitude in the first 
bin as for the full range, but it decreases more strongly 
with centrality and approaches zero for central collisions 
(<5%, (N ch ) — 395). The gap between anisotropies from 
two-pion correlation and conventional flow widens with 
increasing pt as can be seen from Fig.|2Jb). However, the 
statistical accuracy is significantly degraded by invoking 
a two-dimensional window in px- 

The observed excess is attributed to direct pion-pion 
correlations, presumably of semi-hard origin, in addi- 
tion to collective flow. The A<fr distributions shown 
in Fig. 01a,b) are well described by two Gaussians 
at A<j) = 0, 7T on top of elliptic flow. Fit parame- 
ters are the Gaussian amplitudes and widths and back- 
ground B, while V2 is fixed independently by the EP 
method (u|(EP) replacing p in Eq.(2)). The results in 
Fig. 0Ja) show that the close-angle peak stays narrow at 
(To = (0.23±0.03) rad, consistent with fragmentation (lj| . 
The back-to-back peak broadens with centrality up to 
avr = (1.26±0.28) rad at (5-13)%, from where on it can- 
not even be discerned from background. Within the sta- 
tistical errors, the yield contained in both peaks grows 
linearly with N co u which supports the suggested inter- 
pretation of semi-hard scattering. So, the back-to-back 
component escapes detection in central collisions due to 
broadening but does not appear to be suppressed in yield. 

Partonic rescatterings in medium cause an imbal- 
ance i n pt p erpendicular to the initial hard scattering 
plane |19i . l20l ] leading to px broadening that is reflected in 
the width of the back-to-back peak |2l|. In contrast, the 
close-angle peak is not affected since both pions originate 
from fragmentation of the same parton, and propagating 
as color singlets experience only little rescattering there- 
after. In a small-angle approximation, the px broadening 
is (Ap2)V2 a ( PT )(al - cr2)V2. For (5-13)% centrality 
we obtain {Ap T ) 1 / 2 = (1.8±0.4) GeV/c. A more accurate 
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FIG. 5: In-plane (a) and out-of-plane (b) two-pion opening 
angle distributions. Dashed lines are calculated for pure el- 
liptic flow as measured by the EP method and corrected for 
BEC. Data are for centrality 13-26%, pr > 1.2 GeV/c, a cut 
on A9 > 20 mrad, and are efficiency corrected. 



treatment |22j yields (2.8±0.6) GeV/c. Both estimates 
are similar to values measured in pA collisions |2l| . 

Strong non-flow contributions to the two-particle 
opening-angle anisotropy are confirmed by distributions 
in which one of the two pions is detected in the EP 
(±7r/4), or perpendicular to it (again ±7r/4). The 
anisotropy, calculated |2^| for pure elliptic flow as 



Pin/out = V 2 ■ (cos(2(0 - *))) to/out > 



(3) 



is shown by dashed lines in Fig. using the measured 
r>2 =8.5%. The angular brackets indicate an average over 
the respective in-plane and out-of-plane sectors, weighted 
with the anisotropic yields of Eq.(l), and folded with the 
measured EP resolution. 

In-plane (Fig- Ufa)), elliptic flow and semi-hard pairs 
both peak at A<f>= 0, 7r where we observe an excess over 
flow. Out-of-plane (Fig. Eft))), the flow pattern is shifted 
by 7r/2, and the semi-hard correlation (always peaking 
at Acf>= and it) fills in the minima of flow. Would the 
jet-like correlations actually be misidentified elliptic flow, 
then an increased harmonic amplitude in the out-of-plane 
correlation would show more negative swings at and 7r 
than the dashed line in Fig.OJb), which is not the case. 

Summary and discussion.— Differential v 2 (pt) de- 
pends on centrality and px as expected from hydrody- 
namical calculations, but we have not been able to re- 
produce both the magnitude of our v 2 data and the pt 
spectra. Above 1.5 GeV/c the slope of measured v 2 flat- 
tens while the calculation continues to rise. 

The observation of semi-hard two-particle azimuthal 
correlations embedded in collective flow is novel for SPS 
energies. The broadening of the back-to-back correlation 
with increasing centrality suggests that we observe in- 
medium partonic scattering which affects both parts of 
a dijet independently. Although the back-to-back corre- 
lation broadens in central collisions, there is no sign of 



its suppression. The absence of broadening of the close- 
angle correlation supports the view that these pions orig- 
inate from fragmentation of the same parton. Our results 
thus exhibit similar features, but also important differ- 
ences, to recent findings at RHIC |24j . 
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